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ABSTRACT 

We analyze the frequency dependence of the dispersion measure (DM), the column density of free electrons 
to a pulsar, caused by multipath scattering from small scale electron-density fluctuations. The DM is slightly 
different along each propagation path and the transverse spread of paths varies greatly with frequency, yielding 
arrival time perturbations that scale differently than the inverse square of the frequency, the expected depen¬ 
dence for a cold, unmagnetized plasma. We quantify DM and pulse-arrival-time perturbations analytically for 
thin phase screens and extended media and verify the results with simulations of thin screens. The rms differ¬ 
ence between DMs across an octave band near 1.5 GHz ~ 4 x 10“^ pc cm“^ for pulsars at ~ 1 kpc distance. 
Time-of-arrival errors resulting from chromatic DMs are of order a few to hundreds of nanoseconds for pulsars 
with DM ^ 30 pc cm“^ observed across an octave band but increase rapidly to microseconds or larger for 
larger DMs and wider frequency ranges. Erequency-dependent DMs introduce correlated noise into timing 
residuals whose power spectrum is ‘low pass’ in form. The correlation time is of order the geometric mean 
of the refraction times for the highest and lowest radio frequencies used and thus ranges from days to years, 
depending on the pulsar. We discuss the implications for methodologies that use large frequency separations or 
wide bandwidth receivers for timing measurements. Chromatic DMs are partially mitigable by using an addi¬ 
tional chromatic term in arrival time models. Without mitigation, our results provide an additional term in the 
noise model for pulsar timing; they also indicate that in combination with measurement errors from radiometer 
noise, an arbitrary increase in total frequency range (or bandwidth) will yield diminishing benefits and may be 
detrimental to overall timing precision. 

Subject headings: ISM:structure - stars; neutron - pulsars;general - gravitational waves 


1. INTRODUCTION 


Pulsar arrival times with sub- microsecond accuracy ar e required for the detection of nanohertz-frequency gravitational waves 
using p ulsar timing arrays (e.g. Eoster & Backer|1990| l and have payoffs in related area s, such as precision tests of theories of 
gravity (|Kramer et al.|2006|l, determination of neutron star masses (|Antoniadis et al.|20131l, and characterization of microstructure 
in the interstellar medium (^ISM) (|Isaacman & Rankin|1977|l. 


A pulse’s time of arrival (TOA) includes a group delay term, fDM(i^) = ATi/^^DM, where the dispersion measure (DM) is the 
integral of the electron number density along the line of sight, v is the frequency and the constants are AT = and the 

classical electron radius rg. A key step in arrival-time analysis is the removal of the dispersion term. Estimates of DM based on 


TOA differences between two or more frequencies show that DM is epoch dependent for most well-studied pulsars ([Isaacman & 


Rankin|1977l Hamilton et al.|1985[ Cordes et al.|1990 

Phillips & WolszczanI 19911 Backer et al. 19931 K^spi et al. 1994 Ilyasov 

et al.|2005 |Ramachandran et al.|2006||Demorest et al. 

2013||Keith et al.|20131|Petroff et al.|2013 |Eonseca et al. 120141). Temporal 


variations include a slow, systematic trend as a pulsar moves toward or away from us along with stochastic variations from pulsar 
motion that causes the line of sight to sample electron-density fluctuations on a variety of scales. 


In this paper we show that the DM is also frequency dependent because the ISM is sampled differently due to multipath scattering. 
The DM along each ray path is slightly different and the size of the ray-path bundle increases monotonically with decreasing 
frequency, leading to a net difference in dispersive delays in the sum over all ray paths. The variation with frequency of DM(i/, t) 
stays constant for a refract ion time, which is the time scale for the ra y-path bundle to move by its transverse extent and can 
range from hours to years ( Rickett et al.||1984) [Stinebring et al.||2000| l. Our analysis includes the general case of an arbitrary 
medium described by a wavenumber spectrum that varies along the line of sight (LOS). We give specific results for the case of 
a thin screen and for a medium with homogeneous statistics for electron-density variations. We compare analytical results with 
simulations of scattering and dispersion. Our analysis is for the strong scattering regime where all of the flux from a pulsar is 
scattered. In the summary and conclusions section we briefly discuss the weak scattering regime and its role in precision timing. 


Our results complement those of |Lam et al.| ( |2015] l, who quantified TOA errors that result from estimates for DM that make 
use of timing observations at different frequencies that are not made on the same day. Here we show that even simultaneous 
multi-frequency measurements yield TOA errors because the DM varies with frequency. 
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Section 1^ derives the basic effect and gives specific results for a thin screen and for a medium with uniform statistics. Sec¬ 
tion concerns implications for pulsar timing methodology and precision. Sectiondiscusses our results and summarizes our 
conclusions. In Appendix [A| we derive scattering quantities and in Appendix [B] we aerive the RMS two-frequency variation of 
DM. 



Figure 1. Geometries for scattering from a thin screen (left) and filled medium (right). Black lines show simulated ray paths at frequency v and red lines are for 
frequency v' = v 12. The pulsar-ohserver distance is D. The thin screen is at a distance Ds from the pulsar at p and D' = D — Ds from the observer at o. 


2. DISPERSION MEASURE VARIATIONS IN TIME AND FREQUENCY 


There are several underlying causes for temporal variations of DM but only electron density fluctuations in the ISM produce a 
large enough frequency dependence to be important in precision timing]^ For a spatially uniform Ue, the DM is independent 
of frequency and any epoch dependence comes from changes in the pulsar distance D. Motions of the pulsar and Earth with a 


3 X 10 ® pc cm ^Viopt^ 


i- Uen o ver tyr years for a typical 


combined velocity lOOVioo km s ^ parallel to the LOS yield dDM _ 

average electron density Ug = 0.03 cm“^. Linear trends in DM are indeed seen (e.g. Keith et al.^2013 i but this effect is likely to 
be frequency independent. 

Chromatic DMs result from multipath propagation caused by diffraction and refraction from interstellar electron density varia¬ 
tions. The effect we identify is not due to geometrical path-length differences alone. Gravitational lensing, for example, could 
produce multiple ray paths through a medium with constant electron density, but the resulting DM variations would be negligible 
10“^® pc cm“®) given observational bounds on time delays between paths < 1 /is. Moreover, they would be achromatic. 


Density variations exist over a wide range of length scales from kpc to around 10® km (e.g. Armstrong et al.['l995| and the 
smaller ones are responsible for multipath propagation. The DM varies slightly between ray paths and the cross-sectional area 
of the ray-path bundle at any position along the LOS is strongly frequency dependent, cx Density microstructure therefore 
plays two roles: causing multipath propagation and providing variable path integrals of the electron density. Ligure shows 
frequency-dependent ray paths for a thin scattering screen and for a medium that fills the volume between us and a pulsar; it also 
defines some of our notation. 

We model the electron density to be ne(x) = ng{x.) + 6ne{x), where ne(x) is a constant local mean and Sne(x) is the zero-mean 
fluctuating part described by a wavenumber spectrum Psn^ (q, z) that can vary slowly along the line of sight (z-axis). We adopt 
a power-law wavenumber spectrum of the form 


PsuM,z) = Cl{z)q 


qo < q < qi- 


( 1 ) 


The dependence on only the magnitude of q is a simplifying assumption, discussed further in the next subsection, that conflicts 
with some observations that show elliptical scattered images but is consistent with others that show circular images. 

The electromagnetic pha se perturbation cj) from the index of refraction in a cold plasm^is proportional to the integrated electron 
density along a ray path (Rickett|1990[ Appendix A, Eq. A5), 


The effects of the solar wind can be observed in millisecond pulsar observations within 
observations are unsuitable for pulsar timing because of elevated system temperatures and instrumental 
Sun. 


=-Xrg Jdzne{x'{z),z), (2) 

10-solar radii of the Sun jYou et al.|2QQ7} ; however these 
atal distortions associated with the radiation from the nearby 


^ We assume the electron density and magnetic field are small and that frequencies are large enough so that only the term linear in electron density is important. 
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Figure 2. Trajectories of a dispersed pulse in the time-frequency plane for 100 realizations using a phase screen with (pp = 30 rad at 1 GHz. The dedispersed 
trajectories are shown aligned with t = 0, and the inset shows details of the dedispersed trajectories for the lower part of the frequency band. 


where x is the transverse location in the observation plane a distance D from the pulsar while x' is a transverse vector at a location 
z along the LOS. 

For any single ray path the phase perturbation corresponds to a dispersive time delay = {l/2TT){d(j)/dv) characterized by 
a dispersion measure, DM = —^/Arg. Converting Are to DM units at 1 GHz yields a natural value for one radian of phase, 
5DMi = 3.84 X 10“® pc cm“^. Since we know that pulsars are scattered significantly at radio frequencies, phase perturbations 
(j)p on the Fresnel scale rp = ^JXDj^p ~ 10^^ cm are necessarily many radians. The Fresnel scale corresponds to less than 
one day of transverse motion for pulsar velocities ~ 100 km s“^. We therefore expect DM to vary by large multiples of i5DMi 
on time scales of weeks and longer because the RMS phase grows with increasing time span for media with a spectrum like that 
in Eq. By similar reasoning, we expect frequency-dependent DM differences 3> 5DMi because ray paths at widely spaced 
frequencies have separations much larger than the Fresnel scale. 

Figureshows the manifestation of chromatic DMs in simulated arrival times vs. frequency. The curve in the main part of the 
figure IS dominated by the mean DM = 10 pc cm“^. The inset shows the deviations from the mean curve from 0.4 to 0.8 GHz 
for 100 realizations of a phase screen with (j)p = 30 rad, demonstrating the spread in arrival times over ±1 /rs at the lowest 
frequency. 


2.1. Simulations of Phase Screens and Ray-path Averaging 


Phase screens were simulated u sing approaches similar to those presented in Cordes et al. ( |1986| l; Coles et al. ( 1987| l; Foster & 
[Cordes ( |1990| l; [Hu et al.| ( |l991| ) and [Coles et al.|p010| l. An array in the wavenumber domain was hlled with white, Gaussian, 
Hermitian noise shaped by the square root of the ensemble-average power spectrum. A 2D inverse FFT yields the phase screen. 
To include low-frequency components excluded by the FFT, we separately added wavenumber components with periods up to 
three times the array size in each direction. Phase fluctuations were scaled to a specified Fresnel phase at a fiducial frequency, 
which we took to be 1 GHz. We exploited the scale invariance of the phase fluctuations by simulating a set of phase screens for 
small values fp = 5 rad, which require fairly small 2D arrays, and rescaling to larger values. We verified that this procedure was 
correct by explicitly simulating a few large phase screens. 


Figure l^shows an example phase screen and the averaging areas projected onto the screen for a 5:1 range of frequencies. DM 
differences are due solely to the change in averaging area. At different epochs the averaging areas sample different parts of the 
screen, producing stochastic variations in DM that are correlated over the time it takes for the averaging area to move by its 
diameter. The difference in DM between a pair of frequencies is constant over the correlation time for the higher of the two 
frequencies, which has the smaller scattering disk. 

Figure 1^ shows realizations of DM(j^, t) at eight frequencies (left-hand panel) obtained using a phase screen for a Kolmogorov 
spectrum (/3 = 11/3). The figure demonstrates that variations in DM on long time scales are similar over a 5:1 frequency range, 
but differ on short time scales, with a more rapid variation at higher frequencies. The variation in DM has a very ‘red’ power 
spectrum and is a non-stationary process. However, the DM differences between pairs of frequencies (right-hand panel) have 
flatter spectra and show characteristic time scales associated with ray-path smoothing. The autocorrelation widths are roughly 
the geometric mean of the smoothing lengths at each pair of frequencies. These smoothing lengths correspond to the times scales 
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Figure 3. Phase screen for electron-density fluctuations with a Kolmogorov wavenumber spectrum that has 0^ = 5 rad RMS phase variation on the Fresnel 
scale. The x and y coordinates are in units of the Fresnel scale at 1 GHz and the RMS phase difference grows as between two points sepai'ated by 

A = [(Ax)^ -h (Ar/)^]^/^. For a velocity Vx = 100 km s~^, each sample corresponds to 0.28 d and the entire span of x corresponds to 2 yr. Left: phase 
screens without (top) and with spatial smoothing. The bottom three panels are for frequencies of 1.5, 0.75, and 0.3 GHz, as labeled in the right hand panels. 
Plotted circles (barely discernible in the middle two frames) are projected scattering-disk sizes and represent the region with 1/e radius on the phase screen that 
is averaged to produce the DM difference at any given epoch and frequency. The area scales as Right: phase variations con'esponding to trajectories along 

the dashed lines in the left-hand panels. The ID curves are not smoothed versions of (l){x) in the top panel because the smoothing is two dimensional. 


of refractive scintillations. Consequently the time series of DM differences between a pair of frequencies is a ‘red’ process with a 
Gaussian-like power spectrum. Figure [5 shows power spectra of DM variations at two radio frequencies (red curves) and power 
spectra of DM differences for several pairs of frequencies (black curves). 

2.2. Frequency-Dependent Averaging Over Ray Paths 

The measurable DM is an average over all ray paths that reach Earth. At any location along the integration path, the narrow 
ray-path bundle spans a transverse area described by x( 2 ) that also depends on frequency. We define a frequency-dependent 
averaging function Ai^(x, z) that is normalized to unit area, / (ix ^^^(x, z) = 1. The averaged DM is the LOS integral of the 
convolution of the averaging function with the electron density, 

^D nOO 

DM(i^,x) = (DM(x)) + / dz' / dx'Ay(x — x',z')(5ne(x',z'), (3) 

J 0 J — OO 

where the frequency-independent term (DM(x)) is the ensemble average (denoted by angular brackets) of the integrated electron 
density fie- An appropriate averaging function A^(x, z) is the scattered image of the pulsar that generally also includes refraction, 
which offsets the image from the actual pulsar direction. The scattered image 1(0) is the two-dimensional Fourier transform of 
the visibility function r(b) = exp(—D^{b)/2), where is the phase structure function defined in Appendix [ a| The phase 
structure function is the mean-square phase difference of the wavefield between two points separated by distance oTransverse to 
the LOS. In accordance with our assumptions about the spectrum Psn^i D<p depends only on the magnitude of b. 

In the most general case, the image can be elliptical or split into multiple subimages. Elliptical images have been seen for 
some s cattered sources while others are nearly axisymmetric. An extreme case for a nearby pulsar was given by [Brisken et al. 
(|2010[), who infer an image that is Gaussian-like but with ~ 1 % of the flux spread into a long, linear image with high aspect ratio. 




































5 




X (Fresnel Scale Units) 



Spatial Lag (^x (Fresnel Scale Units) 


Figure 4. (Left) A single realization showing DM(t) at 8 frequencies from 0.3 GHz (red) to 1.5 GHz. The DM values have been normalized by the quantity 
for 1/1 = 1 GHz. (Right) Top panel: time series of DM(t) differences between 1.5 GHz and all lower frequencies. The thinnest lines are for the 
frequencies nearest 1.5 GHz and the thickest for 0.3 GHz. Bottom panel: normalized autocorrelation functions of the difference time series in the top panel 
along with the ACF of the DM difference between the two lowest frequencies (0.4 and 0.5 GHz, red). The main lohe of each ACF near zero lag represent the 
characteristic spatial scale of the DM difference. Spatial units can be converted to time units using Ax = Vi At. Using 100 km s“^ velocity and a Fresnel scale 
for 1 GHz and 1 kpc distance, one spatial unit corresponds to 0.14 d. 



Figure 5. Power spectra of dispersion measure variations DM {i/, x) (red lines) and of dispersion-measure differences between frequencies DM(tz') — DM(!z) 
as designated in the legend. The results are based on 1000 realizations of a screen withipp = 30 rad and a Kolmogorov spectrum, which yields a power spectrum 
for DM variations with a —8/3 slope, as indicated by the dashed line. 


Refraction angles are evidently smaller than scattering angles based on the lack of significant angular wandering in VLBI images. 
In addition, for media having wavenumber spectra like those that are consistent with a wide range of scattering and scintillation 
measurements, refraction angles are expected to be small. Since our goal is to describe the basic phenomenon of frequency- 
dependent DMs rather than address all possible varieties of scattering and refraction, the essential features are captured with a 
simplihed approach that allows analytical tractability. The calculation undertaken here simply looks at the difference in electron 
column density that results from the frequency dependence of the ray-path bundle that samples the medium. The calculation of 
the arrival-time difference does not take into account additional delays that result from lateral shifts of the ray-path bundle due to 
refraction or changes due to focusing and defocusing by quadratic phase changes. We emphasize therefore that our results likely 
underestimate the total frequency-dependence of DMs. 

First we make some simple estimates of DM variations. From the phase structure function D^{b) we show that the RMS variation 
of DM across the screen-averaging area scales as the square of the RMS Fresnel phase (pp, and in the next section we derive a 
precise scaling law for the RMS difference in DM between two frequencies. The DM structure function is equal to D^{b)/{Xr^)^. 
Without any ray-path averaging, it can be written in terms of the RMS Fresnel phase as 

4^1 (b\ 

(Are)2 \rF ) 


-Ddm(^) 


(4) 
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Figure 6. Plot of ADM normalized by V(^'^p (u) vs frequency ratio r = v jv' for 100 realizations (black lines). Also shown is -iiFp (r) scaled by the coefficient 
in Eq.|^{red lines), which is the it one standard deviation range of plotted values. 


Averaging over ray-paths yields an RMS DM difference obtained b y averaging th e D M st ructure function over a circular area 


with radius ax 




> rp (which can be derived using Eq. 


A8 


A9 


A12 


andB8 


for moderate to strong scattering with (pp ^ 1. We discuss weak scattering in Election X 


This prescription is applicable only 
The averaging radius is essentially 
the ‘refraction’ scale for refractive scintillations, which is equal to the observed scattering angle projected onto the screen at a 
distance D' from Earth (c.f. Eigure[^; it is also the minimum scale for DM variations. This gives for p in GHz, 


DMr 


(pp/Xr^ ~ 3.84 X 10 “pc cm “ Ppp. 


(5) 


The scaling of DM^ns as the square of (pp arises because the screen phase and the averaging radius are both approximately linear 
in (pp. The Eresnel phase scales with frequency as so DM^ns oc 

In our detailed analysis, we use averaging functions that are symmetric and concentric at different frequencies. In particular, we 
adopt a symmetric Gaussian function for Ai^(x, z), 

Ai,(x,z) = [2TTax{z,p)] ^ ( 6 ) 

where the one-dimensional width ax is proportional to the observed scattering angle and therefore scales with frequency as 
ax oc. with xe ~ 2. We have conhrmed numerically that the Gaussian averaging function yields nearly identical results to 
usage of the scattering image appropriate for a Kolmogorov spectrum of fluctuations. 


2.3. Two-frequency DM Differences 

The difference in DM between two frequencies p and p' (measured at the same location x, corresponding to the same epoch), 

ADM(iz, p' ,yf) = DM(iz', x) — DM(iz, x), (7) 

has an RMS difference, 

/ _ 2 \ 1/2 

= ([^DM(iz,i/', x)] ^ . (8) 

Electron-density wavenumber spectra with (3 > 2 produce DM variations that are dominated by the largest scales. Eor these /3, 
the relevant scales range from the smoothing length ax ~ D'Op^ to the ‘outer scale’ of the spectrum that is likely determined 
by the sizes of structures and clouds in the ISM. The largest structures correspond to time variations in DM up to 10^ years or 
longer for parsec scales and characteristic velocities ^ 100 km s“^. We are concerned with much shorter time scales for which 
the mean-square difference is an appropriate tool because large scale variations ^ ax cancel out. 

In Appendix [Bj we derive the RMS DM difference using Gaussian smoothing functions and express the result in two forms, one 
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that uses the scattering measure SM = f dz C^(z) and a second that uses the Fresnel phase evaluated at the higher frequency, 

f GpQp Te Scattering Measure 




5/3 9/3 


CTe J 


Fresnel Phase. 


(9) 


The dimensionless quantities depend only on the wavenumber spectrum while the dimensionless quantities Gp,gp depend 

on the LOS distribution of Cj^ (thin screen vs. statistically uniform medium, etc.) as well as on the spectrum. Values for 
Kolmogorov media are given in Table All of the relative frequency dependence is contained in the function F’/ 3 (r), 


Ffiir) 


2(4-/3)/2 


^2/3/(/3-2) 


(13-2)12 


-rP - I 


1/2 


( 10 ) 


where r = vjv'. This vanishes for r = 1, as expected and increases monotonically with r for 2 < /? < 4. 
For a Kolmogorov spectrum. 


■ 3.76 X 10-5 pc cm-3 Gn/g £)5/6j,-ii/6 SM.g.g 


^Dm(^:^')=-P'i1/3(3') X 


4.42x10 5 pc cm ^ gwj-i 


v4>%{t2) 


1000 


Scattering Measure 


Fresnel Phase. 


( 11 ) 


where SM is expressed in units of 10 j^pc m fjjg distance D in kpc, and the frequency v in GHz; 4>f{v) is in radians. 

The LOS is characterized by either SM or (j)F, which can be estimated from the scintillation bandwidth and time scale Afiss. 
These are defined respectively as the half-width at half-maximum of the intensity correlation function vs. frequency difference 
and as the half-width at 1/e of the intensity correlation function vs. time lag|^For a thin screen, 

SM_3,5 = 0.74 [(D,/D)(l - D./D)]""/® (12) 

while for a statistically homogenous medium with constant C^, 

SM_3.5 = 2.26(Az/dD)”’^/® . (13) 

To get (/i?’ we use the relationsh ip between the pulse broadening time and the sc intillation bandwidth, = Gi/27rAr'd, where 
Gi is a constant of order unity ( Cordes & Rickett|1998 Lamb^t & Rickett |1999| l, 

-|(/3-2)/4^ ,, X (/3-2)/4 /v/Az/d\^/^^ 


(j)F{l/) = x/2 


Cl 

{D/Ds){l-Ds/D) 




9.6 rad 


100 7 


(14) 


where the quantity in square brackets ~ 0(1) and the approximate expression is for a Kolmogorov spectrum. A ratio v/Avd 
100 is of order the value for nearby millisecond pulsars at z/ = 1 GHz. A similar expression can be written in terms of the 
diffractive scintillation time scale Afjss and the effective velocity 14ff by which the line of sight changes with time. 


~ 8 rad (x/z2l4ff,ioo^iiSS,iooo 


-5/6 


(15) 


with 14ff = 10014ff,ioo km s and Atjss = 10^ Afiss.iooo s. A relation between scattering measure and (j)Fm given in Eq. A9 


Figureshows ADM plotted against frequency ratio r = v jv' (black lines) based on simulations of phase screens with = 
30 rad at 1 GHz. We also show the predicted RMS v') (red lines) indicating statistical consistency between the theoretical 

and simulation results. Individual realizations show a tendency for ADM to be dominated by a linear dependence on r — 1, though 
there are counterexamples that show a significant quadratic dependence. In Sectionj^we utilize this approximate linearity as a 
basis for possible mitigation of the effect in timing data. 


3. IMPLICATIONS FOR TIMING ACCURACY 
There are many complications to the estimation of arrival times. Recent work has addressed the removal of time-variable DMs 


(e.g. Pennucci et al.|20f4l Lee et al.|2014 Liu et al.|2014 1 but no work to our knowledge has aimed at mitigating the frequency 
dependence. Other chromatic ISM effects include delays ffom diffraction and refraction that have been discussed elsewhere (e.g. 


Foster & 


Cordes|1990||Rickett|l 

1 pulsar’s pulse (e.g.lCr; 


19901 ICordes & Shannon|2010|l. Also important is the evolution with frequency of the intrinsic 

“ TtfigTOl 


shape of a pulsar’s pulse (e.g.|Craft & Comella|1968(|Craft|l976l|Ahuja et al.|2007||Hassall et al.|2012[|Pennucci et al.|2014||Liu 

et al. 2014|l, which produces a systematic TOA error vs. frequency that is partially covariant with DM v^iations but is assumed to 


® Usage of the half widths at half maximum and at 1/e for the two cases is natural given the mathematical forms of the correlation functions for media with a 
square-law phase structure function. 
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be independent of epoch. Since these and other effects are largely independent, their variances add and therefore can be discussed 
separately from the frequency-dependent effect analyzed here. 


We consider multiple frequency measurements at a fixed epoch, so we drop any explicit time dependence from DM. A simple 
model for the TOA at frequency v includes dispersive delays and measurement errors as perturbations of the ‘true’ TOA, 


a:dm(j^) 


(16) 


The quantity is an additive, frequency-dependent error that is uncorrelated between different frequencies for radiometer-noise 
but is highly correlated for intrinsic pulse jitter, at least over modest frequency separations. 


3.1. Dual Frequency Measurements 


An observing strategy exploits the lever-arm provided by widely-spaced frequencies to obtain a high-precision estimate for DM, 
which then is used to estimate a DM-compensated arrival time, too - Here we derive the resulting errors in DM and too. 

Consider two arrival times measured at frequencies v and v' < v the same epoch. The usual operational practice is to 

estimate DM (denoted by the caret) by inverting Eq. [funder the assumption that DM is frequency independent, 

tu ’ tu 


DM = 

K(v'-^ - ^-2) 

The dispersion delay is then removed from the measured TOA to estimate the infinite-frequency TOA, 

^ ^ KJm _ ^ ^ K[UM{v) - i 5 m ] 

too - tiy n - loo T o T 

V V 

For a frequency ratio r = vjv', the difference between the estimated and true DM(z/) is 


DM(z/) - DM = 


[DM(i/) - DM(i^')] + 


and the resulting error in the infinite-frequency TOA is 

_ K 

atoo — loo loo - o 


-I 


) ABM{v,v') + 


K 


e^' - r-e,, 

j-2 _ 1 


(17) 


(18) 


(19) 


( 20 ) 


The estimator too is unbiased if the DM variations and the errors e,^ have zero mean values over an ensemble. The contribution 
to 6too from the measurement error at the higher frequency is enhanced by a factor r^. This means that the error at the lower 
frequency / can be larger than the high-frequency error by a factor equal to a modest fraction of and not unduly affect the 
precision of too ■ 

The combined RMS timing error cry is the quadratic sum of the individual RMS errors. 




1/2 


( 21 ) 


3.1.1. TOA Error from Frequency-dependent DMs 
The RMS DM difference defined in Eq.[9|T0| implies an RMS error in TOA, 

7-2 Ka 


{2tt) SM, Scattering Measure 




|r2-l| 


T- -=Epir) X 


5/3 7/3 




27ric 


( 22 ) 


Fresnel Phase, 


where the r-dependent factors have been combined into a timing-error function valid for 2 < /3 < 4 (see Figure [T3]l, 


E/sir) = 


|r2-l| 


( 23 ) 


For a Kolmorovov spectrum, the RMS timing error is 

{ 156 ns Gii /3 SM_ 3,5 Scattering Measure 

184 ns 7 ii /3 711/3 




1000 




Fresnel Phase. 
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Figure 7. RMS error in the DM-corrected TOA too vs. frequency for a fixed upper frequency of 2 G Hz and four values of average DM, as labeled . Th e curves 
for different D Ms were calculated using the pulse-broadening time from the scaling law in Eq. 7 of|Bhat et al.H2004| to obtain (ppi^) using Eq.|14|and then 
evaluating Eq.[^for different r = vjv'. The spread in values demarked by the shaded regions results trom tne variation about the empirical mean relation 
between log;i^Q and log^Q DM. While we show curves extending down to 0.1 GHz, pulse broadening will likely dominate the timing precision and will render 
some pulsars undetectable at lower frequencies. 




Frequency (GHz) Frequency (GHz) 


Figure 8. RMS too vs. frequency v' for dual-frequency observations for pulsars with DM = 10 pc cm~^ (left) and DM = 30 pc cm“^ (right). For each curve 
the rightmost frequency is the highest frequency of the pair and the curve gives crt^ (ly, v'). Lower v have higher crtoo (^5 ^ ')■ Low, mid, and high-frequency 
bands are indicated. 


Figure 0 shows sT)M{u' v = 2 GHz for four values of average DM. We calculated these curves by estimating 

the Fresnel phase at the reference frequency v = 2 GHz using Eq. and then expressing the scintillation bandwidth at the 
highest frequency v in terms of the pulse broadening time using = Ci/27rTrf. We get Td from the empirical relation 


Bhat et al. 


2004 1 . The variation about this 


logio '''difJ-s) = —3.46 + 0.154x + 1.07a:^ — [2/3/(/3 — 2)] log^o for v in GHz (e.g. 
mean relation is 0.7 in log^Q Tdins). We use /3 = 11/3 though Bhat et al. find a best-fit value that is slightly greateiQ The 
curves in the figure may therefore underestimate the timing error because we use v = 2 GHz. The figure demonstrates the strong 
dependence of the TOA error on chromatic DMs because Td and thus (j>F increase nonlinearly with DM and the TOA error scales 
as the square of 


Figure [^complements Figure |7| by showing the TOA error for observations that involve different frequency ranges and for two 
different but rather low values of DM (10 and 30 pc cm“^). Each curve in the figure corresponds to two-frequency observations 
where the higher frequency v is the right-most frequency (2, 0.8, and 0.4 GHz) and the lower frequency v' is the horizontal 
axis. Obtaining TOAs with < 0.1 ps precision clearly requires a high frequency of no less than 0.8 GHz for low-DM pulsars. 


^ Specifically, Bhat et al. fo und a best fit a = 3.86 it 0.16 for the coefficient of the logj^Q u term, which corresponds mathematically Xa fd = la j(a — 1 = 
4.15 it 0.17, but this expression applies only for /3 < 4. So a conservative interpretation is that Bhat et al.’s result corresponds to ^5 > 11/3. Alternative 
interpretations involve other aspects of the scattering geometry rather than the form of the wavenumber spectrum. 
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The curves presented here and in Figure [^demonstrate that low frequency telescopes such as LOFABj^need to be coupled with 
higher-frequency telescopes operating at ^ 2 GHz to provide TOA errors ^ 0.1 ps from chromatic DMs. Note however, that 
the total DM-corrected timing error also involves radiometer noise and pulse jitter, which favor wide frequency spans that are 
discussed next. 


3.1.2. Measurement Errors from Radiometer Noise, Pulse Jitter, and Diffractive Scintillations 

To treat any kind of measurement error, we use a correlation function (normalized to unit maximum) 
between the errors at two frequencies. The RMS of the second term in Eq.j^is then 

[r^cr^+ cr2, _ 

-• (25) 


Radiometer noise (e = rn) is uncorrelated between TOAs obtained using non-overlapping bandpasses, so ' = 0 and 


s/r^al + al, 

|r2-l| 


(26) 


In our analysis we consider wide frequency ranges with logarithmic spacings of frequencies and assume that bandwidths are 
proportional to frequency. This naturally leads to a power-law scaling of the TOA error from radiometer noise. 


We then rewrite Eq.[2^as 


CTrn.iy' 


O-r-n.i/ (r'^ + r 


> 1/2 


|r2-l| 


(27) 


(28) 


Individual pulses show phase and amplitude jitter (e = j) that cause small shape changes in the averages of large numbers of 
pulses used to calculate TOAs. The resulting TOA error is correlated between frequencies, sometimes highly so. If perfectly 
correlated {p^^i = 1) and identical, jitter produces no error in DM because the TOAs move in tandem at the two frequencies 
(c.f Eq. [T^. The resulting TOA error from jitter is then simply Uj. However, single pulses and average profiles evolve slowly 
with frequency, yielding random and systematic TOA errors, respectively. Generally the jitter correlation will be less than 100%, 
yielding a larger contribution to the TOA error than from perfectly correlated jitter. However, the frequency s eparation over which 
single pulses decorrelate is large for the few cases that have been studied. These include the Crab pulsar ( [Sallmen et al.|l999) 
which decorrelates ov er about 0.7 GHz and the millisecond pulsar J0437-4715, which decorrelates over a 2-GH z bandwidth 
( Shannon et al.||2014|l. Eor the brightest pulsars, the TOA errors from jitter and radiometer noise are comparable ( [Shannon &| 
|Cordes|2012[|E)olch et al.|2014| [Shannon et al.|20f4] l so when pulse jitter is not completely correlated, mis-estimation of DM is 
comparable to that from radiomieter noise. 


We adopt a power-law scaling analogous to that for radiometer noise. 




= r 


The resulting expression for the jitter-induced timing error for perfect correlation between frequencies {pi, 

o -i.u I 


= 1) is 


r2-l 


(29) 


(30) 


Diffractive scintillations (DISS) have correlation times and bandwidths that range, respectively, from minutes to hours and ^ 
100 kHz to 100s of MHz at 1 GHz for currently timed millisecond pulsars (MSPs), which tend to have low DMs ^ 30 pc cm“^. 
DISS causes TOA errors because the associated pulse broadening function — the scattering impulse response of the ISM that 
is convolved with a pulsar’s pulse — is stochastic on the DISS correlation time scale. The RMS TOA error at an individual 
frequency is much smaller than those from radiometer noise and jitter for nearby MSPs but can be comparable for more distant 
ones (Lam et. al, in preparation). 

Scintillations of a low-DM MSP will yield a non-zero correlation p^^/ for observations made nearly simultaneously (e.g. within 
less than one hour) and with frequencies separated by no more than a correlation bandwidth. Most current timing observations 
are made with large-enough bandwidths or frequency separations and many observations are made non-simultaneously, in which 
case p^v' = 0. 

* http://www.astron.nl/radio-observatory/astronomers/technical-information/lofar-technical-information 
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r r 

Figure 9. (Left) RMS error of too for two-frequency measurements as a function of r = v' jv for a fiducial (highest) frequency of 1 GHz. Individual 
contributions are shown from radiometer noise, pulse jitter, and frequency-dependent DMs along with the total (heavy black line). For this case we assume that 
the noise and jitter errors are both 100 ns and do not depend on frequency. Results are not dramatically different for cases where they are frequency dependent. 
(Right) Total TOA error vs r for six different mixtures of timing errors from radiometer noise, jitter, and frequency dependent DMs. The legend gives the TOA 
errors from RMS noise crt,rn, from RMS jitter, and their frequency scaling indices, Xm and Xj respectively. The dashed curve gives the TOA error if only 
radiometer noise with 10 ns RMS contributed to timing errors. 


3.1.3. Systematic Errors from Profile Evolution 


Profile evolution (e = ’pe’) is known to introduce systematic errors in TOAs because the average profile changes smoothly with 
frequency. They can also be described using Eq. l20| However, unlike the random errors in the previous section, profile evolution 
can be mitigated by exploiting its apparent epocnindependence (e.g. Pennucci et al.|20l4) l and thus differs from the chromatic 
DM effect that varies with epoch. 


3.1.4. Assessment of Two-Erequency Timing 


It is often assumed that wider frequency separations yield more precise dispersion measures and arrival times. We assess this 
approach by considering the timing error from frequency-dependent DMs combined with measurement errors. The results indi¬ 
cate that there can be diminishing if not worsening returns once the frequency ratio becomes larger than about a factor of two 
and if no mitigation procedure is used. We note that the same affliction arises from evolution of pulse shapes with frequency, as 
mentioned earlier. 


In Figure |9j we plot af from Eq. |2T and the individual contributions to it from frequency-dependent DMs, radiometer noise, 
and pulse jitter using Eq. 24 2^ am 30 respectively. The specific case shown uses RMS values of 100 ns for both the noise 
and jitter contributions that areTrequency independent (aim = Xj =0). The DM error was calculated for a scattering measure 
logj^o = —3.5, a value that is typical of pulsars within 1 kpc of Earth. The figure shows basic features that are generic. Rather 
than decreasing monotonically or remaining flat, as expected from radiometer noise or pulse jitter alone, the TOA error reaches 
a minimum at r ^ 2 and then rises as the frequency-dependent DM contribution begins to dominate. Other cases with different 
mixtures of radiometer noise and jitter and different indices Xm and xj 0 are qualitatively similar. 


Other examples are shown in Figure]^ including those with very small (10 ns) contributions from noise and jitter. While some 
cases — those with large noise or jitter contributions — do not show a minimum in r, all show much larger asymptotic TOA 
errors than expected from radiometer noise and jitter alone. This implies that increases in bandwidth with wideband systems will 
provide diminishing returns unless the frequency-dependent DM can be mitigated, as discussed below. 


3.2. Wideband Timing Measurements 


Digital backend systems developed recently for pulsar observations have much larger frequency ranges (1.8:1) than previous 
systems and provide the opportunity to e stimate TOAs over a c ontinuous range of frequencies rather than using two narrowband 
frequencies with a wide separation (e.g. Pennucci et al.|[20T4 l. Even wider bandwidth systems are contemplated with 4:1 (or 
larger) frequency ranges. 


We analyze wideband systems with arbitrary total bandwidths by using standard least-squares fitting methods to fit data without 
and with frequency-dependent DMs. We let the data vector D comprise a set of TOAs k = 1,..., N,j} and their mea¬ 
surement errors for N^, separate frequencies measured at the same epoch. For simplicity, we consider only radiometer noise 
in the wideband analysis, corresonding to a diagonal covariance matrix, C = diagonal For a design matrix X and a 
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Figure 10. RMS TOA error from wideband fitting for DM and too over a continuous frequency range from 0.4 to 2 GHz. The phase screen has a Fresnel 
RMS phase of 30 rad and the RMS radiometer noise is an optimistic 10 ns (left) and 100 ns (right) that is assumed to be constant in frequency. A total of 2000 
phase-screen realizations was used. The cases shown include (1) Fixed true DM is fixed and a timing model that includes only too and a term; (2) Variable 
DM without noise added and a fit for fixed DM; (3) Variable DM with noise and a fit for fixed DM; (4) As with (3) but with a fit that also includes a term 
with fixed JV = 3; and (5) As with (4) but where the exponent X is also fitted for. 
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Figure 11. Results from wideband fitting for DM and too over a continuous frequency range from 0.4 to 2 GHz. The phase screen has a Fresnel RMS phase of 
30 rad and the RMS radiometer noise is constant in frequency: 10 ns (left set of panels) and 100 ns (right). Left column: histograms of the error in too based on 
2000 realizations of a phase screen. RMS values for too are given in each panel. Right column: histograms of the difference in DM from the value at the highest 
frequency. RMS values for (5DM are given in units of 10”“^ pc cm“^. First row: the true DM is fixed and the timing model includes only a term; Second 

row: same as the first case except that DM varies with frequency; Third row: DM varies with frequency and the timing model includes both a and a 
term with JV = 3 (fixed); and Fourth row: same as previous except that the exponent is also fitted for. 


linear model D = X0 + e, the solution vector is 0 = (X^^C ^X) ^ X^^C and the covariance matrix for the parameters is 

Vq = (XtC-iX)”\ 

With respect to the timing model of Eq. we consider four alternatives; 


1. First is where the frequency-dependent DM is negligible so DM is constant in frequency and the only errors are from 

radiometer noise. This would be the case for a very low-DM pulsar or measurements at high frequencies 1 / ^ 1 GHz. The 
data are htted with a parameter vector 0 = col (too, ^t'DM) and corresponding design matrix X = matrix(l i = 

1,..., The solution vector is unbiased. 

2. The second case is where the frequency dependence of DM is signihcant but the data are still htted with a constant DM 
model. The results are generally biased. 

3. The third case is motivated by recognizing in Figurej^that DM differences ADM(i/', u) have a tendency to appear roughly 
(but not exactly) linear in r — 1; this suggests that a term in the htting function scaling as will absorb much of the 
effect and improve the estimate for too- 
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4. The fourth case is the same as the third except the fitting function includes ■&. v ^ term where x is also fitted for as a 
(nonlinear) parameter instead of fixing it at x = 3. 

Example results in Figure show the TOA error plotted against total bandwidth for :/ = 2 GHz (the highest frequency) with 
a total of 20 separate frequencies used. DM variations are for a phase screen with = 30 rad at = 1 GHz. In the left- 
hand panel we have used an optimistic 10 ns for the radiometer noise at the fiducial frequency of 1 GHz while in the right-hand 
panel we have used 100 ns. We do not include pulse jitter in these examples because results shown previously indicate that it is 
secondary. Curves are shown for the four cases described above along with a fifth case that is an extreme version of case 2 where 
the noise is assumed negligible. The results are similar to those obtained in the two-frequency case. Namely, for the 10-ns noise 
case, increases in bandwidth yield improvements until the bandwidth is about 1 GHz and then the results degrade if there is no 
explicit fitting for the frequency-dependent DM. With such fitting (cases 3 and 4 above), the bandwidth can be increased another 
20 to 40% before the timing errors degrade. For the larger 100-ns noise in the right-hand panel, TOA errors are dominated by 
radiometer noise except for the largest bandwidths for which there is an increase in timing error. 

Figure shows additional information on fitting results using histograms of the errors in too and DM. The two cases shown 
are for tHe same fits and data model used to produce Figure [T^ The top rows show fitting results when the DM is constant 
in frequency and the sole source of fitting error is the additive noise. Other rows show the results for the different methods 
outlined above to deal with the frequency dependence of DM. Overall, Figures and [TT] show that allowance for the frequency 
dependence can reduce the timing error but cannot achieve the same results as tot a constant DM. 


3.3. Other Mitigation Methods 


So far we have discussed estimation and removal of dispersion using data obtained at a single epoch. To date, most methods used 
by groups aiming to detect gravitational waves have used multiple epochs t o estimate DM at the epoch of any particular arrival 
time.While using multiple epochs can be deleterious (e.g. Lam et al.|2015| l, they can be implemented with algorithms that take 
into account the correlation times of DM variations that result from ray-path averaging. These times can be days to months or 
longer depending on the frequency and ISM along the line of sight. 


It is beyond the scope of this paper to develop a multi-epoch approach. However, we can illustrate a two-frequency approach 
that smooths the DM time series at a high frequency to match the characteristic time scale of the low-frequency DM variations. 
Figure [T^ shows DM variations at 0.2 and 1.5 GHz from a simulation for a phase screen with (j)p = 30 rad at 1 GHz. In the 
top panek the high-frequency time series has been optimally smoothed using a Gaussian smoothing function that minimizes the 
mean-square difference with the low-frequency time series. The bottom panel of the figures shows the smoothed and unsmoothed 
DM difference between the two frequencies. While the differences have been reduced by smoothing the high-frequency DM, 
they are not negligible because one-dimensional smoothing cannot model the two-dimensional smoothing that occurs in the ISM 
from scattering. Nonetheless, this approach can be used as a mitigation procedure. 
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Figure 12. Test of one-dimensional smoothing to estimate the two-dimensional averaging from scattering using DM variations for a scattering screen with 
= 30 rad. The left-hand axis label is DM units and the right-hand axis is in time-delay units. (Top): DM(z/) vs. x for two widely space frequencies, 0.2 
and 1.5 GHz. Also shown (red line) is the smoothed 1.5 GHz variation that best matches the 0.2 GHz curve using a Gaussian smoothing function. (Bottom): DM 
differences without and with smoothing of the high frequency variation. 
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3.4. Weak Scattering 


Our analysis applies to the strong scattering regime where (l)p > 1. Nearby pulsars observed at higher frequencies may be in weak 
scattering with (/ip- < 1. In this regime, the scattered pulsar image consists of the very compact image of the unscattered pulsar 
that contains a fraction ~ 1 — of the total flux combined with a scattered image with the remaining fraction. For small (f>F, 
the cross-sectional area that is averaged is very small and the DM will be nearly achromatic. This would be another advantage of 
using high frequencies — defined as th ose a bove the transition frequency to weak scattering. The transition frequency is derived 

8 GHz For some pulsars, the transition frequency is 


by requiring (pp < 1 rad and using Eq. 


A9 


to get Vtr 


^ 1 GHz (e.g. Rickett et al.|2000 1 . However, precise DM estimates require multiple-frequency (or wideband) observations that 
almost always will use lower frequencies that are in the strong-scattering regime. We therefore see no way to avoid the chromatic 
aspect of DMs in precision timing. 


4. SUMMARY AND CONCLUSIONS 


We have shown that dispersion measures are chromatic because microstructure in the interstellar electron density causes multipath 
propagation that is strongly frequency dependent. We have characterized the effect in terms of the average over ray paths, 
DM, using an averaging kernel that is frequency dependent. Results were given for media having a power-law electron-density 
wavenumber spectrum and for arbitrary variations in amplitude of the spectrum along the line of sight. We verified our analytical 
results with simulations of phase screens with a Kolmogorov spectrum. Differences in DM between two frequencies scale as 
oc SM, where pp is the RMS phase across a Fresnel scale at the highest frequency used and SM is the scattering 
measure. Observations at frequencies ^ 1 GHz will have pp ^ 1 for most pulsars and ipp oc The specific variation of 

DM with frequency will remain constant in time over a refraction time scale of hours to months (or longer) that is generally larger 
for larger-DM pulsars. The joint variation of DM in time and frequency thus differs from that of intrinsic profile evolution, which 
appears to be epoch independent for millisecond pulsars. Longer period pulsars show state switching (nulling, mode changes, 
etc.) on a wide range of time scales, but the switching statistics appear to have stationary statistics. 


As yet, there has been no direct detection of chromatic D Ms. This is not surprising be cause there has been no detailed analysis 
prior to the present paper other than a brief discussion in |Cordes & Shannon] pOlO)!, and its effects would be impossible to 
identify uniquely in dual-frequency TOA measurements typically employed. However, departures from the dispersion delay 
appropriate for a tenuous, cold, unmagnetized plasma have been searched for and use d to put limits on chromatic timing effects 
that scale differently than the standard dispersion law (e.g. Hassall et al.||2012 and references therein). Such departures 
are complicated by profile evolution with frequency that is mitigated by identifying flducial pulse phases that yield consistency 
with the cold-plasma law. Most of the objects analyzed this way (pulsars and fast radio bursts of apparent extragalactic origin) 
have much coarser timing precision than the MSPs used in pulsar timing arrays. The best prospects for detection are from a 
bright, high-DM pulsar that has minimal profile evolution over the frequency range needed to probe DM variations. Profile 
evolution may be disentangled from chromatic DMs by exploiting the lack of (or minimal) epoch dependence of profile evolution 
in comparison with variations in DM that have a characteristic correlation time for each frequen cy. It is conceivable that some of 
the frequency-dependent timing variations observed from the MSP J1909—3744 (e.g. Fig. 7 in [Manchester et al.|2013| l include 
the chromatic DM effect. 


Chromatic DMs have a significant impact on pulsar timing applications where sub-microsecond timing accuracy is needed, 
such as detection of gravitational waves with pulsar-timing arrays and high-order general relativistic effects in binary pulsars. All 
pulsar timing applications depend on removing dispersion delays with high accuracy to estimate TOAs unaffected by propagation 
through the ISM, which we have called too- We have analyzed errors in too resulting from methodologies that assume DMs 
are achromatic for two cases, one where TOAs are measured at pairs of widely separated frequencies and a second that uses 
continuous, wideband systems. Chromatic DMs introduce errors in too that depend strongly on the pulsar’s mean DM and on the 
particular range of frequencies used. For nearby pulsars with DM ^ 30 pc cm“^ and an octave frequency range with 1.5 GHz 
as the highest frequency, TOA errors solely from the chromatic part of DM are of order a few to hundreds of nanoseconds for 
observations extending down to 1 GHz or 0.2 GHz. Timing errors increase rapidly with increasing mean DM and timing residuals 
will be correlated on time scales related to those of refractive interstellar scintillations. 


Timing errors from radiometer noise and pulse jitter combined with chromatic DMs show a broad minimum as a function of 
total frequency range (or bandwidth) that is of order an octave in frequency. This arises because TOAs improve monotonically 
with bandwidth as far as radiometer noise is concerned, but the opposite is true for chromatic DMs. The simplest prescription 
for optimizing TOA precision is to use an upper frequency that is as high as possible. The choice of upper frequency is strongly 
pulsar and telescope dependent. 


Chromatic DMs also need to be considered in combination with other chromatic effects, including intrinsic pulse profile evolution 
with frequency and additional interstellar delays that result directly from scattering and refracti on that have different frequ ency 
dependences than any of the DM effects. A comprehensive assessment of effects like that injCordes & Shannon] (|201()|) that 
includes the chromatic DM effect is deferred to a separate paper. Even though the timing errors from chromatic DMs are smaller 
than other effects, they nonetheless may inhibit improvements in timing accuracy that otherwise might be obtainable. It is 
therefore important to confirm that chromatic DMs are present in timing data at predicted levels and develop ways to mitigate 
them, if possible. If not, they need to be part of the noise model for timing analysis and incorporated into the covariance matrix 
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used in model fitting. 

In a separate article we will assess the role of chromatic DMs for all objects currently being observed in pulsar timing array 
campaigns to detect gravitational waves and we will also assess different methodologies for using existing and future telescopes. 
We will also identify particular pulsars that are good candidates for direct detection of DM chromaticity. 
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work reported here. JMC’s work on pulsar timing at Cornell University is supported in part by NSF PIRE program award number 
0968296. RMS acknowledges travel support through a John Philip early career research award from the CSIRO. DRS gratefully 
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APPENDIX 


A. PHASE STRUCTURE EUNCTION AND SCATTERING ANGLE 


The phase structure function is the integral from a point source to an observer at distance D from a point source, 

r ^ r 

D4,{h) = (^(t){yL) - (t){yL + b)f'^ = AT:{r^Xf j dz y dqj. (qj., z) (^1 - 

where b is a spatial offset and qj^ a vector wavenumber transverse to the line of sight. This form applies to a wavenumber 
spectrum Psn^ whose extent in wavenumber is much narrower than D~^, has a shape independent of z, and an amplitude that 
varies slowly with z. Normalization is so that the mean-square electron density is the integral of Psn over all wavenumbers, and 
C^dzP^^^^ = 27rS{q.). 

Scattering measurements indicate various degrees of anisotropy of density fluctuations, but the isotropic case is easier to analyze. 
To treat the anisotropic case is tedious and does not add any further insights to the results obtained for the isotropic case. For 
isotropic irregularities only the magnitudes of b and qj^ matter, yielding 


= 87r^(reA)^ / dz 

Jo 

where Jq is the Bessel function of the hrst kind. We adopt a power-law wavenumber spectrum, 

PsnA<ij.,z) = Cl{z)qr^, qo<q<qi, (A3) 

where C^(z) varies slowly with z on length scales much larger than the outer scale, 27r/go- For b (i.e. b smaller than 

the inner scale), the phase structure function is quadratic in b while for 5 ^ it asymptotes to twice the total variance of the 
phase. In the intermediate regime where <C 6 <C gjj"^, 

D^b) « //3(Are)"SMeff b'^-^ (A4) 

where (Cordes & Rickett 1998; Eq. B6) 



J dqj_qj_PsnAdJ-,z)[^-Joid^bz/D)], 


(A2) 


U = 


Stt^ r(2-/3/2) 


(^-2)2/5-2 r(/3/2) ’ 

where the effective scattering measure is the line-of-sight weighted integral of C^, 


SMeff= / dzc: 




^ ^ = SM X = SM ((z/D)^-2) 


/ dz C2(z) 


C5 


(A5) 


(A6) 


where the angular brackets denote an average over the LOS using C^(z) as a weighting function and the scattering measure is 


.D 

SM= / dzC2(z). (A7) 

Jo 

For a screen, SMes = (Ds/D)^“^ SM and for a uniform medium with = constant, SMes = SM/(/3 — 1). A plane wave 
incident on a foreground scattering medium so that Dg/D —^ 1 gives SM^ff = SM (e.g. for an extragalactic pulse incident on 
the Milky Way). Values for SMes/SM are given in Tabie[^for a Kolmogorov spectrum along with other parameters quantities. 

Alternatively we can express the phase structure function in terms of the RMS phase (j)p across a Fresnel radius rp in the screen. 
We dehne the Fresnel scale using rp = [(AiA)/27r]ZJeff/J9 where for a thin screen = DgD'/D = D{Ds/D){1 — Dg/D). 




16 


For a statistically uniform medium we assume Dgs = D/A. Taking into account that a transverse scale b at an observer’s position 
corresponds to a scale zb /D at a position z along the LOS, we have for a thin screen at z = Dg, 




Dg 

D 


From this and Eq. A7 we solve for SM in terms of c/p, 

SM = 


(A8) 


(A9) 


(Are)^' 

In the following we assume the same relation holds generally though we have derived it from the thin-screen case. 

The scattered image of a p oint source has longer tails than a Gaussian function for /3 < 4 and an inner scale ^p/qi much smaller 
than the Fresnel scale (e.g. Rickett|1990| l. However it is convenient to characterize the main part of the image with an equivalent 

)n has the 


Gaussian whose visibility function 


fie same 1/e width. This defines the spatial scale bg using D^{bg) = 2, 


b,= 


ll/(/3-2) 


from which the RMS angular size ae, the 1/e half width, and the full width at half maximum (FWHM) are calculated as 


ae = 


ZFWHM 


A 


(AlO) 


(All) 


2 V2 In 2 V‘2TTbe ' 

The power-law spectrum yields an RMS angular size that we factor into a scattering size ag^ and a geometry-dependent factor. 


erg = erg 


^ SM y 


= - 
TT 


A^Z/jr^SM 


2/5/2 


l/(/3-2) 


(A12) 


B. TWO-FREQUENCY CROSS CORRELATION AND STRUCTURE EUNCTION 


We calculate the mean-square of the difference ADM(i^, i/', x) = DM(i^', x) — DM(z^, x) defined in the main text to get the 
two-frequency structure function, 

= ( [ADM(i^, I/', x)]^ = Cg^{F, v) -f - 2C^^{v', v) (Bl) 

where the cross correlation of i5DM = DM — (DM) between two frequencies is 

= (i5DM(i/,x)5DM(i^',x)) = JJ dx'dx." JJ dz'dz" A^(x', (x", z")(5ne(x', z')fce(x", z")).(B2) 

The z integrals are from 0 to D and the x integrals are over an infinite plane. We define the cross-correlation function of the 
averaging function. 


Ca{Sx, = j dx — (5x/2, z)Ay(a: — (5x/2, z), 


(B3) 


and assume that it changes slowly in z. By changing variables from x',x" to x = (x' -f x")/2 and (5x = x' — x" and from 
zz" to z = (z' -f z")/2 and 5z = z' — z", and using the hierarchy of scales assumed in Appendix [ a] the integration over 6z 
gives 2Tr5{qz) and we obtain 

= \ J dzj dSyi J dSzDsn,{d^,dz;z)[2CAiS:!^,z,v,F') - CAiS:>e,z,v,v) - CAiS:!^,z,v',ly')]; (B4) 

Dsn^ ((5x, 6z; z) = ([5ne(x, z) — 6ne{x + (5x, z -I- Sz)]'^'/ is the structure function for the electron density. For the power-law 
spectrum of Appendix [A| we have 


J dSz Dsn^ (Sx, Sz, z) = ffiCliz) (Sxf ^ . 


(B5) 


We now adopt the Gaussian smoothing function of Eq.l^and employ the frequency scaling of cr^, which is the same as that of the 
RMS scattering angle erg oc with xe = /3/(/3 — 2j. By referencing frequency-dependent quantities to the higher frequency 
F we get 




( 2 ^) 2 r (2 - /?/ 2 ) ^2 


/3-2 


J dzCl{z)a^ ^(z,f) 


(B6) 
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Figure 13. (Left) The frequency scaling function -F^(r) vs. r = vjv' for three values of /3. (Right) The frequency scaling function E^{r) for the same three 
values of 0. 


where the frequency-scaling function is 


i 2(4-/3)/2 ri+^2/3/(/3-2)' 


(/3-2)/2 


-rP - I 


1/2 


(B7) 


Figure [o shows Fp{r) along with a related function Ep{r) = r‘^Fis{r)/(r‘^ — 1) that is used in the main text to characterize 


time-ot-airival errors. 


To evaluate the integral in Eq. B6 we need expressions for how the transverse extent of the ray path bundle varies with z. We 
define 


H. = 


ax{z,iy) = Dae{v)h{z/D) 
1 — x) for a statistica’ 

J dzCl(z)[h{z/D)f-^ _ 


(B8) 

where h{x) = 1 — x for a thin screen, h{x) = x{l — x) for a statistically homogeneous medium, and h{x) = 1 for a plane wave 
incident on the medium. We then define 


/ dzCliz) 

istically unil 

incidence, Flp is the limit of the thin-screen case as Dg —> oo. 


^{[h{z/D)f-^) 


Cl 


(B9) 


For a thin screen, Flp = {D'/D)d for a statistically uniform medium, F[p = dw [r(;(l — w)]^ and for plane-wave 


Using Eq. A12 to evaluate ae{v), the RMS DM difference becomes 

^dm(^> = G/30/3 77(/5-2)/2 


where the two dimensionless quantities are 


Q/3 = 


{V2TT)^-dT{2-(3l2)fp 


1/2 


Gb = 


Hi- 


/SMeff\ 

V SM y 


-I 1/2 


(/3-2) 

All of the geometry-dependent factors are consolidated into Gp. 

We also express v') in terms of the Fresnel phase by substituting for SM from Eq. A9 


(BIO) 


(BID 




7/3 = 


2^/27r2r(2 -/3/2) 


1/2 


90 = 


Hb 


Values of G/j, qp and are given in Table[T]for a Kolmogorov spectrum. 


/ D /SMeff\ 
\d,s) [ SM J 


1/2 


(B12) 
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Table 1 

Relevant Factors for a Kolmogorov Medium (/3 = 11/3) 


Geometry Independent Faetors: 



Quantity 

/3 = 11/3 



fp 

88.3 



Qp 

22.0 



<ip 

1.15 



Fp(^) 

■^ 11 / 3 ( 1 ) = 0 

Fn/3(2) = 1.056 

Fii/3(5) = 6.64 

Geometry Dependent Factors: 



Quantity 

Thin Screen 

(x = Ds/D) 

Uniform 

Plane Wave 

DefflD 

x{l — x) 

1/4 

1 

SMeff/SM 

3,5/3 

3/8 

1 

h(z) 

1 - z/D 

(z/D){l-z/D) 

1 

Hp 

(1 — 

0.056 

(1 — x)®/® as X ^ 0 

Gp 

[x(l-x)]®/6 

0.145 

1 

9P 

1 

0.46 

1 
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